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Abstract 

The RAP experiment was recently carried out in order to explain the 
results in cosmic ray detection obtained by the gravitational wave resonant 
antenna NAUTILUS when operated above and below the temperature of 
transition from superconducting to normal state of the material. In this 
experiment, performed by shooting an electron beam on a suspended bar 
made of Al 5056 alloy (~ 5% Mg), both the amplitude of the bar first longi- 
tudinal mode of oscillation, excited by the beam interacting with the bulk, 
and the energy deposited by the beam in the bar were measured. These 
quantities, inserted in the equations describing the mechanism of the mode 
excitation and complemented by an independent measurement of the specific 
heat, allow us to determine the linear expansion coefficient a of the mate- 
rial. We obtain a = [(10.9 ± 0.4)T + (1.3 ± 0.1)T 3 ] x 10" 10 K" 1 for the 
normal state of conduction in the temperature interval 0.9 < T < 2 K and 
a = [(-2.45 ± 0.60) + (-10.68 ± 1.24)T + (0.13 ± 0.01)T 3 ] x 10~ 9 K+ 1 for 
the superconducting state in the interval 0.3 < T < 0.8 K. 



"Corresponding author. Tel.: +39 06 94032622; Fax: +39 06 94032256. Email Address: 
carlo.ligi@lnf.infn.it 



Preprint submitted to Physica B 



Keywords: Thermal expansion, Low temperature, Aluminum alloy 
PACS: 65.60.+a, 61.82.Bg, 41.75.Fr 



1. Introduction 

Wrought aluminum-magnesium alloys (International Alloy Designation 
System: 5XXX, Mg content: 5.6% maximum) are commonly used in ap- 
plications where good workability, very good resistance to corrosion, high 
fatigue strength are desidered. Example of applications are: oil, fuel lines 
and tanks, pressure cryogenic vessels, marine structures and fittings, auto- 
motive and architectural components. 

The alloy Al 5056 (5.2% Mg, 0.1% Mn, 0.1% Cr) is the material of 
the test mass of gravitational wave (GW) resonant antennas AURIGA [l[, 
EXPLORER Q, NAUTILUS Q in Italy and ALLEGRO 3 in the United 
States, operated at liquid Helium temperatures and below. During the last 
decade, all these detectors took part, as a network, in coordinated searches 
for impulsive GW excitations (no GW events were detected) jsl, 0]. 

The test mass of the GW resonant antennas is a suspended cylinder 
(~ 2300 kg mass), acting as a mechanical resonator whose resonances are 
excited by the incoming GW. The elastic vibrations are converted to elec- 
trical signals by a transducer system. The minimum detectable energy of a 
GW resonant antenna inversely depends on Q, the acoustic quality factor of 
the material also known as the inverse of the internal friction. Al 5056 was 
chosen as the material for the GW antennas due to the high values of Q 0] , 
exhibited at very low temperatures. 

B.L. Baron and R. Hofstadter first measured mechanical oscillations 
in piezoelectric disks penetrated by high energy electron beams and they 
outlined the possibility that cosmic ray (CR) events could excite mechanical 
vibrations in a metallic cylinder at its resonant frequency and that these 
events could mimic GW events detected by resonant antennas. NAUTILUS 
has been equipped with a CR detector [i[ to study the interactions caused 
by CRs and to provide a veto against the CR induced events in the antenna. 
Correlating NAUTILUS data with CR observations, the following results 
were obtained: 1) the rate of high energy signals due to CR showers was 



larger than expectations [lOj, [ll| when the antenna was operated at a tem- 



perature T = 0.14 K, i.e. well below the transition temperature from normal- 



conducting (n) to superconducting (s) states of Al 5056 12]; 2) this feature 
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was not observed when the antenna was operated at T = 1.5 K, i.e. above 
the transition temperature 13 . 

Since these findings had not a straightforward interpretation, we performed 
an experiment (RAP) aimed at measuring the longitudinal oscillations of a 
small suspended Al 5056 cylinder hit by an electron beam of known energy. 
The RAP experiment, performed at the Beam Test Facility (BTF) [l4| of the 
DAFNE ^-factory complex in the INFN Frascati Laboratory, has shown that 
1) the amplitude of the longitudinal oscillations of the cylinder hit byioniz- 



ing particles depends on the state of conduction of the material [15|, [16J and 
that 2) the observed amplitudes are consistent with the amplitudes measured 
by NAUTILUS in CR events both in in) and (s) states [171 ] . Moreover, the 
amplitude can be evaluated in the framework of the Thermo- Acoustic model: 
in the (n) state it depends on the ratio of the thermal expansion coefficient 
to the specific heat, this ratio being part of the definition of the Griineisen 
parameter. In the (s) state the amplitude also depends on the fractional 
volume change between the (s) and (n) states. In principle, the relevant 
thermophysical properties of the Al 5056 at low temperature are not known 
with sufficient precision to allow the comparison between experimental and 
theoretical values of the amplitude predicted by the model. However, our 
confidence in the model is also based on the small discrepancy, of the order 
of 15%, between the amplitudes measured in the alloy and those computed 



by using the data for pure aluminum in the region 4.5 <T< 264 K [18 



Direct measurements of the Griineisen parameter of elastic materials have 



been reported [19|: the authors used pulsed electron beams hitting the front 
surface of a thin slab, and measured the induced motion on the rear face sur- 
face. The equation of motion describing thermoelastic effects depends on the 
Griineisen parameter and on the energy deposited by the beam in the slab. 
The measurement of the displacement, the knowledge of the deposited energy 
and the use of the equation of motion allowed them to directly determine the 
Griineisen parameter. 

In a similar way this paper presents experimental results on the linear ex- 
pansion coefficient of the aluminum alloy above and below the temperature 
of transition between the (s) and (n) state. In order to determine the coef- 
ficient we use: a) the oscillation amplitude measurements, b) the equations 
describing the amplitude in the framework of the Thermo-Acoustic model, 
c) the measurements of the deposited energy by the electron beam, obtained 
by the product of the measured beam multiplicity and the calculated energy 
loss in the material by each particle, and d) an independent measurement of 
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the specific heat. 



2. The Thermo- Acoustic model 

2.1. (n) state 

We consider a cylinder (radius R, length L, mass M) that is suspended 
in correspondence of the middle section, with its axis horizontal. A ionizing 
particle, after hitting the lateral surface of the cylinder and interacting with 
the material, generates a pressure pulse in the bulk. This sonic pulse is the 
result of the local thermal expansion of the material caused by the warming 
up due to the energy lost by the particle in the material. The sonic pulse de- 
termines the excitation of the vibrational oscillation modes of the suspended 
bar. By introducing a vector field u(x, t) to describe the local displacements 
from equilibrium, we express the amplitude of the /c-th longitudinal mode of 



oscillation as proportional to the quantity |20|, [21|, 122 



^Ptherm ai-j- 
9k,therm ^--L-k 
P 

dW 



dl 



(1) 



where AP t herm is the pressure pulse due to the sonic source described above, p 
is the mass density, dWj dl is the energy loss per unit length of the interacting 
particle, A' is the cross section of the tubular zone centered on the particle 
path in which the effects are generated and = J dl(V ■ Ufc(x)) is a line 
integral over the particle path involving the normal mode of oscillation Ufc(x). 
In the previous relation 7 is the adimensional Griineisen parameter: 

pcv 

with (3 the volume thermal expansion coefficient ((3 = 3a for cubic elements), 
cy the isochoric specific heat and Kt the isothermal bulk module. 

A thin bar (R/L 1) hit by a ionizing particle at the center of the lateral 
surface is a particular case of the more general one that leads to relation (CQ) . 
In this case, the fundamental longitudinal mode of oscillation is excited to a 



maximum amplitude given by 23 



x - 2aLW (3) 
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and then exponentially decays due to internal friction. Here W is the total 
energy released in the bar by the particle beam. The experimental conditions 
of RAP are close to this particular case and the amplitude maximum of the 
fundamental longitudinal mode of oscillation is the observable measured by 
the experiment. However, in order to model the amplitude in the most 
realistic way, we have performed a Monte Carlo (MC) simulation [l8j, which 
takes into account the corrections 0[(R/L) 2 ] for the modes of oscillation of a 
cylinder, the transverse dimension of the electron beam at the impact point 
and the trajectories of the secondary particles generated by the electron 
interactions in the bar. All these effects are summarized in a corrective 
parameter e and the amplitude maximum of the fundamental longitudinal 
mode of oscillation is modeled by: 

X (n) = X{1 + e) (4) 

The value of e for the aluminum alloy bar used in the experiment is estimated 
by MC to be -0.04. 

2.2. (s) state 

The energy released by the ionizing particles to the material determines 
the suppression of the superconductivity in a region (Hot Spot) centered 
around the particle path. The maximum value of the Hot Spot radius, rns, 
is obtained by equating the energy lost per unit length by the particle, dW/dl, 
to the enthalpy variation in the volume undergoing the (s) — > (n) transition 
at temperature T (2^, 25, 26]: 




\dW/dl\ 

rns = \ jr-r 1 = \ — , (5) 

V 71 All V 7T 

where Aft, is the enthalpy change per unit volume and A" is the cross section 
of the zone switched to the (n) state. 

The amplitude of the longitudinal oscillations is the sum of two terms, one 
related to thermodynamic effects in the (s) — > (n) transition and the other to 
the thermal effects in the region switched to the (n) state, which have been 
already described in the Section 12.11 The contribution to the amplitude 
of the cylinder k — th oscillation mode due to the (s) — > (n) transition is 



proportional to [20|, |21 



A P 

t—i± trans a/It 
Qkjrans ^ -*-k 
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where V = V( s ) and AV, AS are the differences of the volume and the 
entropy in the two states of conduction. According to the thermodynamics 
of volume and pressure effects for Type-II superconductors, these differences 
can be written as 00: 



AV = V (n) - V (s) 



and 



AS 



«S(ra) — <S( S ) 



VH C dH c 
4tt dP 



VH c dH c 
4vr dT 



HldV 



HldV_ 
8^df 



(6) 



(7) 



Here H c is the superconducting critical magnetic field that is supposed to 
have the parabolic behavior H c (t) = H C (0)(1 — t 2 ), where t = T/T c and 
T c is the transition temperature. Moreover, the terms in { } brackets are 
expected to be smaller than the preceding terms jiif and are usually ignored 
for practical purposes. 

Finally, the knowledge of the specific heat of the material for the (s) state, 
c^ s \ allows us to compute the variation of the entalpy per unit volume in the 
form 291: 



Ah = C int (T) + T 



AS 



(8) 



with: 



Ci n t{T) 



(T')dT' 



The amplitude maximum of the fundamental longitudinal mode of os- 
cillation in the (s) state, as already indicated, is the sum of two terms, 
one related to the transition and the other to the (n) contribution from the 
switched region: X (s) = X trans + X (n) . By observing that X trans /X( n ) = 
g k=o, trans/ 9 k=o, therm we can describe the amplitude maximum according to 
the following relation: 



1 We show the relations in cgs unit system, as done by the original authors. In the 
computations we will convert the magnetic field energy density from cgs to SI. 
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V V 

where the relation and the definition (J2J) of 7 are used to define: 

2pL{\ + e) 



n 



3irM(X (n) /W) 



We note that the term contained in { } brackets in the relation (Q is inde- 
pendent from W, the energy released by the particle to the bar, and that 
X( s ) linearly depends on W through X( n y 

3. Experimental setup and procedures 

The experiment mechanical layout, its cryogenic setup and operations, the 
electron beam characteristics, the instrumentation and the procedures for cal- 



ibrations, data taking and analysis have been fully described in Refs. [171 . |18 
Briefly, the dimensions and the mass of the Al 5056 cylindrical bar are 
R = 0.091 m, L — 0.500 m, M = 34.1 kg, respectively. The bar hangs 
from the cryostat top by means of a multi-stage suspension system providing 
a 150 dB attenuation of the external mechanical noise in the 1700-6500 Hz 
frequency window. The frequency of the fundamental longitudinal mode of 
oscillation of the bar is fo = 5413.6 Hz below T = 4 K. The cryostat is 
equipped with a dilution refrigerator. The temperatures are measured inside 
the cryostat by 11 thermometers controlled by two multi-channel resistance 
bridges and, among them, a calibrated RuC>2 resistor measures the temper- 
ature of one of the bar end faces with an accuracy of 0.01 K below T = 4 K. 
Two piezoelectric ceramics (Pz), electrically connected in parallel, are in- 
serted in a slot cut in the position opposite to the bar suspension point and 
are squeezed when the bar shrinks. In this Pz arrangement the strain mea- 
sured at the bar center is proportional to the displacement of the bar end 
faces. The Pz output is first amplified and then sampled at 100 kHz by an 
ADC embedded in a VME system dedicated to the data acquisition. A Pz 
calibration procedure, performed before each run of data taking, provides 
the factor converting the ADC samples into the displacements of the bar end 
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faces. A software filtering algorithm, known as " digital lock- in" , extracts the 
Fourier component at the frequency /o from the time sequence formed by 
the ADC samples before and after the beam hit, determining the amplitude 
of the induced fundamental oscillation. The sign of the amplitude is taken 
positive or negative according to the sign of the first sample raising above 
the threshold after the beam hit. 

The BTF beam line delivers to the bar single pulses of ~ 10 ns duration, 
containing N e electrons of 510 ±2 MeV energy. N e ranges from about 5 x 10 7 
to 10 9 and is measured with an accuracy of ~ 3% (for N e > 5 x 10 s ) by 
an integrating current transformer placed close to the beam line exit point. 
The MC simulation, discussed in Section [2J estimates an average energy lost 
(AE) ± a AE = 195.2 ± 70.6 MeV for a 512 MeV electron interacting in 
the bar and, consequently, the total energy loss per beam pulse is given by 
W = N e (AE), aw = v^Al &ae- Two sources of error affect the vibration 
maximum amplitude: the first is an overall systematic error of the order of 
±6%, that accounts for the slightly different set-up and calibration proce- 
dures implemented in the runs over two years and the second is related to 
the noise in the measurement of the oscillation amplitude (±1.3 x 10~ 13 m). 



4. Linear expansion coefficient measurements 

The coefficient a is determined by inserting the measured values of the 
amplitude maximum of the fundamental longitudinal mode of oscillation 
(X( n ), Xr s )) and the corresponding measurements of the deposited energy 
(W) into the relations modeling the amplitudes (equations (@J,fl9])). The 



measurement [30] at very low temperatures of the specific heat of an Al 5056 
sample, of the same production batch of the RAP bar, is a fundamental 
ingredient for the coefficient determination both in (n) and (s) states. More- 
over, the superconducting characterization of the material has shown that 



T c = 0.845 ± 0.002 K with a total transition width of about 0.1 K [30jj. Thus 
we determine a for the two states of conduction («(„), «( s )) according to this 
temperature value. 

4-1. (n) state 

Values of the amplitude maximum of the fundamental longitudinal mode 
of oscillation and the related values of the deposited energies by the beam 
have been measured and analyzed in the temperature interval 0.9 < T < 
2 K [17]. The values of Xr n ) show a strict linear correlation with the values 
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of W (Fig. [T]) and they have been fitted, including the errors, by a line 
constrained to the origin (Xt n \ = poW). The slope of the fit, which has a x 2 
per degree of freedom equal to 0.45/15, is po = (2.24 ± 0.05) x 10~ 10 m J -1 . 
From the relation @ it follows that: 

o_ = nM Po 

c v 2L(l + e) v ; 

In the same temperature interval the parametrization cy = CT + DT 3 , where 
C is the electronic specific heat coefficient per unit volume and D is the lattice 
contribution, is determined by C = (4.382 ± 0.117) x 10" 2 J kg" 1 K~ 2 
and D = (5.20 ± 0.37) x 10~ 3 J kg" 1 K -4 . Inserting this cy parametrization 
in relation fIXOjl directly gives: 



a 



in) = «(n),e + «(n),<? 



[(10.9 ± 0.4)T + (1.3 ± 0.1)T 3 ] x KT 1U K 



—i 



(11) 



for the linear expansion coefficient (Fig. |2]) expressed in terms of the electronic 
and lattice components in the temperature interval 0.9 < T < 2 K. 



We note that the value of ai( n ) for pure aluminum is reported [31] to be 
2.7 x 10~ 9 K" 1 at T = 2 K. This value, which is the one associated to 
the lowest temperature in the previously cited work, can be compared with 
a( n ) = (3.2 ± 0.2) x 10~ 9 K _1 obtained in the present work for Al 5056 at 
the same temperature. 

4-2. (s) state 

From eq. (Q we can derive the volume variation among the two states of 
conduction: 



AV 
V 



1 

if 



n Ah 

Po 



T 



AS 

V 



(12) 



where 1Z = X^/W — p , with p introduced in Section I4~T1 The volume 
variation depends on measured quantities (X^, W, p , cy, T, T c ) and calcu- 
lated ones. The latter are related to the critical field H C {T = 0) and to its 
parabolic dependance on the temperature, which also allows us to compute 
the derivative dH c /dT. If the superconducting properties of Al 5056 can 
be described by the BCS theory, then H C {T = 0) w 2.42V^T C rj 70 Oe, C 
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being the electronic specific heat coefficient per unit volume (Section 14. ip . 
here expressed in suitable cgs units (erg cm -3 Kr 2 ). 

As already mentioned, the model predicts that the amplitude maximum of 
the fundamental longitudinal mode of oscillation for both states of conduction 
(X(n), Xf a -\) linearly depends on the energy deposited by the e~ beam in 



the bar. On the contrary, the Xr s \ values for Al 5056 |15|, [17J show an 



increasing deviation from linearity with the increase of the deposited energy 



W, due to saturation effects discussed in Ref. [17] . However, if we restrict 
our analysis to the data gathered at the lowest deposited energies, these 
non linearities can safely be neglected, and we can use the relation ( IT2|) to 
compute the volume variation. Fig. [3] shows the values of AV/V versus the 
temperature T for the data with released energy: 1.5 xl0~ 3 < W < 9xl0~ 3 J. 
AV/V data have been fitted by a 2nd order polynomial go + Qi T + q 2 T 2 
obtaining q = (-1.737 ± 0.054) x 10~ 8 , q x = (7.334 ± 1.785) x 10" 9 K _1 , 
q2 = (1.766 ± 0.187) x 10~ 8 K~ 2 with a x 2 — 1-38 normalized over 94 degrees 
of freedom. 

The use of q and of the relation (ED gives dH c /dP = (-3.12 ± 0.10) x 
ig-12 rp Pa -1 at T = K, a value that can be compared with (—2.67 ± 
0.06) x 10~ 12 T Pa -1 obtained by Harris and Mapother in their experimental 
study of the critical field of pure aluminum as a function of pressure and 
temperature 
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The difference of the thermal expansion coefficients in the normal and 
superconducting phases is obtained by taking the derivative with respect to 
T of the difference of volumes |27f : Aa = ai( n ) — a( s ) = f ^jf ^r- 



Under the hypothesis that is still represented by the relation ( ITT]) in the 
interval 0.3 K < T < T c and by using the polynomial fitting for AV/V, the 
following expression is obtained for that temperature interval: 

= «(n) - Aa 
= [(-2.45 ± 0.60) + (-10.68 ± 1.24)T 

+ (0.13 ± 0.01)T 3 ] x 10~ 9 K" 1 (13) 

This result can be compared with the expansion coefficient of pure aluminum 
in the (s) state, using for pure aluminum the relation o(( s ) = Q( s ), e + ate, in 
which it is assumed that the lattice component does not depend on the 
conduction state and that its values are given in Ref. [33] . The definition ([2]) 
gives a( s ), e = pcv, e ^l(s),e/ (3-Kt) > where the electronic component of cy in 
the (s) state (cy^) is taken from the work of Phillips [34], the electronic 
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component of the Griineisen parameter in the (s) state (7( s v e ) from Ref. 35 

The comparison 



and K, 



K, 



(n) 



79.4 x 10 Pa at low temperatures [36 



is shown in Fig. H]as a function of the reduced temperature t. 



5. Conclusions 

In this paper we present an evaluation of the linear expansion coefficient 
of an Al-Mg alloy obtained by measuring both the amplitude of the funda- 
mental mode of the longitudinal oscillation excited by electrons interacting 
in a suspended bar and the energy released in the bar by the electron pulse. 
The use of this method to determine a is new, although similarities exist 
with the direct measurement of the the Griineisen parameter performed by 
hitting thin slabs with particle beams to generate thermoelastic pulses. No 
expansion data for aluminum alloys were previously available in the litera- 
ture for the temperature range 0.3 < T < 2 K explored by this experiment. 
The expansion coefficient is negative in the superconducting state and its ab- 
solute value just below T c is larger by an order of magnitude than the value 
above this temperature. This feature is expected also for pure aluminum, 



according to an analysis |35| of the measured values of H C (P,T) available in 



Ref. 32 
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Figure 1: Normal state. Values of the amplitude maximum of the fundamental longitudinal 
mode of oscillation (X( n \) vs. the deposited energies (W) measured in the temperature 
interval 0.9 < T < 2 K. Data have been fitted by a straight line, constrained to zero, with 
slope (2.24 ±0.05) x 10~ 10 m J" 1 . 
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Figure 3: Superconducting state. Measured values of AV/V vs. temperature for beam hits 
deposing in the bar an energy W < 9 mj. The thick dotted line is a fit to the data with a 
2nd order polynomial. The thinner lines enclose a region compatible with the uncertainties 
in the fit parameters. 
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Figure 4: Superconducting state. Values of ctr s ) vs. reduced temperature (t = T/T c ) for 
Al 5056 (solid line). For reference, we also show the result for pure Al (dotted line) 
derived from data of Ref. fffjj /. The thinner lines enclose the regions compatible with the 
uncertainties in the ai s -\ values. 
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